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^ ' 1 INTRODUCTION 
•rH . 

■ Pulsar-like stars keep to manifest surprising observational 
$_i features since their first discovery in 1967. One of their ex- 
5^ i traordinary behaviors, whose origin and astrophysical impli- 
cations are hotly debated about, is the superflare from soft 
7-ray repeaters (SGRs). Only three such events have been 
detected in three of the four SGRs: 1979/03/05 of SGR 0525- 
66 (spin period P = 8.1 s), 1998/08/27 of SRG 1900+14 
(P = 5.16 s), 2004/12/27 of SGR 1806-20 (P = 7.45 s). The 
released energy in the former two flares could be ~ ^q 44 ~ 4 5 
erg, while ~ 10 47 erg in the third. A peculiar nature of 
superflare is the initial brief (~ 0.2 s) spikes of 7— rays 
with energies up to several MeV, which contain most of 
the fl are energy and are followed by tails lasting minutes 
(e.g.. Ulurlev et all BOPS'). Additionally, quasi-periodic os- 
cillations (QPOs) during superflares were found soon af- 
ter the onsets of the superflares in SG R 1806-20 (at fre- 
quency / ~ 93 Hz. Ilsrael et alJ|2005[) . in SGR 1900+14 
(/ ~ 84 Hz, IStrohmaver fc Wattj|200ai . and in SGR 0526- 
66 (/ ~ 43 Hz. iBarat et al.lll983Fl . Higher frequency oscilla- 
tions at about 150, 625, and 1,8 40 Hz are also detected from 
the superflare of SGR 1806-20 (IStrohmaver fc WattsllioOrj . 

Current models for the superflares are in the scenario 
of magnetars, a kind of neutron stars with polar magnetic 



ABSTRACT 

Three times of supergiant flares from soft 7-ray repeatres are observed, with typical 
released energy of ~ io 44 ~ 47 er g_ A conventional model (i.e., the magnetar model) for 
such events is catastrophic magnetism-powered instability through magnetohydrody- 
namic process, in which a significant part of short-hard 7-ray bursts could also be the 
results of magnetars. Based on various observational features (e.g., precession, glitch, 
thermal photon emission) and the underlying theory of strong interaction (quantum 
chromodynamics, QCD), it could not be ruled out yet that pulsar- like stars might 
be actually solid quark stars. Strain energy develops during a solid star's life, and 
starquakes could occur when stellar stresses reach a critical value, with huge energy 
released. An alternative model for supergiant flares of soft 7-ray repeatres is presented, 
in which energy release during a star quake of solid quark stars is calculated. Numer- 
ical results for spherically asymmetric solid stars show that the released gravitational 
energy during a giant quake could be as high as 10 48 erg if the tangential pressure is 
slightly higher than the radial one. Difficulties in magnetar models may be overcome 
if AXPs/SGRs are accreting solid quark stars with mass ~ (1 — 2)M Q . 
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fields in the range 10 14 ~ 10 15 G, including the SGRs and 
the anomalous X-ray pulsars (AXPs). The quiescent X-ray 
emission with luminosity ~ io 34 ~ 36 erg/s as well as the su- 
perflares of SG Rs are supposed to be powe red by magnetic 
field decay (e.g.. lWoods fc Thompsonl2005lh However, there 
are still some debates on magentars (see §3 for more discus- 
sion on their existence), although they are really popular 
in the astrophysical society. Alternatively, we propose here 
that a giant quake in a solid quark star may result in a su- 
perflares and reproduce the general observational features. 

A quark star is composed dominantly by quark mat- 
ter which is a direct consequence of the asymptotic freedom 
nature in quantum chromodynamics (QCD) , the underlying 
theory believed for the elementary strong interaction. It is 
worth noting that this stellar quark matter at low tempera- 
ture should be very different from the hot quark matter to 
be searched in relativistic heavy ion colliders. A degenerate 
Fermi gas of quarks is expected at extremely high density 
and temperature, while Cooper pairing of quarks near the 
Fermi surface occurs in cool and dense (but not asymptot- 
ically) quark matter because of the strong and attractive 
QCD quark-quark interaction. A condensate of the pairs 
may then result in a color supercond uctivity (CSC) phase in 
this case (lAlford fc Raiagopal2006l and references therein). 
However, a solid state with quark-clusters in periodic lat- 
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tices was also conjectured in a parametric region where the 
density could be lower than that of CSC state 1x112003. 
2005). This hypothetic state could still not be ruled out by 
simple QCD principles as well as astrophysical observations 
iXull200 6a. and references therein). Unfortunately, due to 
the non- linear nature of QCD, one can not now obtain with 
certainty the critical parameters (e.g., baryon density and 
temperature) for the asymptotically free quark phase, the 
CSC phase, and the (solid) quark-cluster phase. 

Pulsar-like stars, including SGRs and AXPs, are com- 
pact remnants of evolved stars, the nature of which is 
still a matter of controversy. Though these stars are pop- 
ularly thought to be normal neutron stars, no convincing 
work, either theoretical from first principles or observa- 
tional, has excluded the possibility that they are actually 
quark stars. Besides the rotational, thermal, magnetic, and 
accretion-gravitational ones, the remaining free energy for 
a solid quark star in cludes additional elastic and quake- 
gravitational energies l)Horvathll200l IxlfeoOfi'l. In this pa- 
per, we will show that these additional free energies of solid 
quark stars could be high enough to power the superflares of 
SGRs. Though the energy budget quest of the superflares is 
focused on, the radiative mechanisms as well as astrophysi- 
cal implications of such events are discussed too. 



2 THE MODEL 



dient, |dP/dr|, decreases if P± > P, which may result in 
a higher maximum mass of compact sta rs. Whether the 
2.1Mo-millisecond pulsar iNice et al.ll2005fl in a binary sys- 
tem with a helium white dwarf secondary is relevant to this 
nature of solid stars could be an interesting topic. One can 
also see that a sudden decrease of P± in a star may cause 
substantial energy release, since the star's radius decreases 
and the absolute gravitational energy increases. 

A quark star would initially be in a fluid state well ap- 
proximated by a perfect fluid. In a simplified version of the 
bag model for quark matter, the equation of state for fluid 
quark matter could be 



(2) 



where the bag constant B ■<? could be between 60 MeV/fm 3 
and 110 MeV/fm 3 . Since no equation of sate for solid quark 
matter is available, we may just approximate the radial pres- 
sure by Eq.((2j in the following calculations. Also the equa- 
tion below holds, 



4irr p. 



(3) 



dm(r) 
dr 

We have to integrate numerically these differential equations 
of Eq. lH4L which are complete for P, p and m(r) if B and 
e are certain, in order to know the e-dependent global struc- 
ture (e.g., radius and mass) of a solid quark star. 

Substituting Eq. @ into Eq.@ and Eq. Q, we have 



It is very fundamental to study static and spherically sym- 
metric gravitational sources in general relativity, especially 
for the interior solutions. The Tolman-Oppenheimer-Volkoff 
(TOV) solution is only for perfect fluid. However, for solid 
quark stars, since the local press could be anisotropic in 
elastic matter, the radial pressure gradient could be partially 
balanced by the tangential shear force although a general un- 
derstanding of relativis tic, elastic bodies has unfortun ately 
not been achieved fe.g.. lKarlovini fc Samuel sson 2004). The 
origin of this local anisotropic force in solid quark stars could 
be from the development of elastic energy as a star (i) spins 
down (its ellipticity decreases) and (ii) cools (it may shrink) . 
Release of the elastic as well as the gravitational energies 
would be not negligible, and may have significant astrophys- 
ical implications. 

Let's numerically calculate the structure of solid quark 
stars as following. For the sack of simplicity, we only deal 
with spherically symmetric sources in order to make sense of 
possible astrophysical consequence of solid quark stars. By 
introducing respectively radial and tangential pressures, P 
and Pj_ , the stellar equilibrium eq uation of static anisotropic 
matter in Newtonian gravity is llHerrera fc Santobs!ll997t 
Eq. (2.4) there): dP/dr = -Gm(r)p/r 2 +2(P ± -P)/r, where 
p, m(r) and G denote, respectively, mass density, mass in- 
terior the radius r and the gravitational constant. However, 
in Ei nstein's g ravity, this equilibrium equation is modified 
fe.g.. lLfi]IT999|) . 



dP = Gm(r)p( 1 + ^)( 1 +#)?) | 2e p 

y2 ^ 2Gm(r) y 



(i) 



where P± = (1 + e)P is introduced. In case of isotropic 
pressure, e = 0, Eq.Q turns out to be the TOV equa- 
tion. It is evident from Eq.Q that the radial pressure gra- 



d r 
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dp Gm(r)(4S+ i |f ) <1+ (4B + 

dr 
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The precision of numerical solution should be very high in 
order to obtain the small-e-dependent features from Eq.lQJ. 
Initially, the core of a star is supposed to be homogeneous, 
with the density p = po, the radius r = 0.1 cm, m(0.1) = 
4po7rr 3 /3, and P(0.1) = (po — 4P)c 2 /3 in the computation. 
Numerical method can then be used to integrate Eq.JlJ from 
r = 0.1 cm to the boundary, r — R, of the star, step by step 
(with index n). At the star's boundary, P(R) = and M = 
m(R), with R = r„ if and only if P{R n ) ^ 0, P(P„+i) < 0. 
Approximated in a flat spacetime, the total gravitational 
energy E and the stellar moment of inertia I could also 
been obtained in the numerical process, by integrating the 
following equations, 



dE 



Gro(r) 



dm(r) 



-^W 2 pdr) 



3P 

= -47r(— + 4B)Gm(r)rdr, 
cr 



" 2 

dl — —r dm(r) 



-r 2 (47rr 2 pdr) 



8tt , 3P , 4 , 
( +4B)r 4 dr. 



(5) 



(6) 



The Runge-Kutta method of order-4 is used in the code. 

Some issues are worth noting during the numerical pro- 
cess. It is very necessary to improve the precision in order 
to obtain the differences, which we are interested, of E(e) 
and 7(e) by numerical integration, since both the related 
numbers and the domains in this problem vary in 7 orders 
of magnitudes. Computing error comes mainly from the de- 
cision of a star's boundary. To improve the precision of in- 
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Figure 1. The difference of stellar radius as a function of stellar 
mass. Solid lines are for bag constant B = 60 MeV/fm 3 , and 
dashed lines for B = 110 MeV/fm 3 . The quantity of e is defined 
by introducing P± = (1 + e)P, where Pj_, P are the radial and 
tangential pressures, respectively. 
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Figure 2. Same as in Fig.l, but for the difference of gravitational 
energy. 



tegration, we have to divide the domain to be smaller and 
smaller near the boundary, until the quantities calculated 
are credible. 

Our calculation results are shown in FigQ to Fig|3] 

Starquakes may result in a sudden change of e, with an 
energy release of the gravitational energy as well as the 
tangential strain energy. Generally, it is evident that the 
differences of radius, gravitational energy, and moment of 
inertia increase proportionally to stellar mass and the pa- 
rameter e. This means that an event should be more im- 
portant for a bigger change of e in a quark star with higher 
mass. It is shown in Fig^that the stellar radius varies in- 
significantly when small e is considered (-R(e) — R(0) is only 
from 10~ cm ~ 10 cm for e = 10" 8 — 4 if stellar mass 
is ~ Mq). Only gravitational energy release is considered 
in Fig|3 which could be as high as 10 48 erg if e ~ 10 -4 
and M ~ Mq. Typical energy of ifj 44 ~ 47 erg is released 
during superflares of SGRs, and we may then propose a gi- 
ant starquake with e < 10~ 4 could produce such a flare. A 
sudden change of e can also result in a jump of spin fre- 
quency, Afi/fi = —AI/I. We may expect from Fig[3]tliat 



Figure 3. Same as in Fig.l, but for the ratio difference of moment 
of inertia. 



glitches with Afl/fi ~ 10 10 ~ 4 could occur for parameters 
of M = (0.1 ~ 1.4)M Q and e = 10" 9 ~~ 4 . It is suggestive 
that a giantflare may accompany a high-amplitude glitch. 



3 CONCLUSIONS AND DISCUSSIONS 

We suggest alternatively that the superflares of soft j- 
ray repeaters could be the results of giant quakes of solid 
quark stars. Numerical calculations for spherically asymmet- 
ric solid stars show that the released gravitational energy 
during a giant quake could be as high as 10 48 erg if the tan- 
gential pressure is slightly higher (only ~ (1 + 10 -5 ) times) 
than the radial one in a star with mass ~ Mq. However, a 
detail process by which the released energy is transformed 
into radiation is still not clear during a quark. The quake- 
calculation presen ted in this paper is different from that of 
IZhou et al.l j2004) . 

A direct consequence in the model here is that, be- 
sides giant quakes, much small quakes especially of aged 
quark stars (both their magnetospheric and thermal radia- 
tions would be low enough not to be detected regularly by 
recent facilities) with low masses are expected to occur as 
hard X-ray transients. There is a rough tendency that both 
the jump amplitudes and the frequency of the glitches are 
apparently decreased with the pulse period when the pulsars 
become old llLvne et alj EoOO'l. No glitch pulsar with the pe- 
riod longer than 0.7 second has been detected. According to 
the law of seismology, no big quake occurs if many small ones 
happen frequently, but a giant quake may take place after 
a long time of silence (i.e., no quake period). If this law ap- 
plies, an old solid quark star with long spin period (S> 1 s), 
which may have enough time to leave its host galaxy, could 
have a big quake (crash) after a long time of silence. Part 
of short 7-ray burst could probably rare flares of middle or 
low amplitudes of starquakes in the galactic halo. More such 
events could be recorded by the Swift or the future HXMT 
(hard X-ray modulation telescope). 

Anomalous X-ray pulsars/Soft 7-ray repeaters 
(AXPs/SGRs) are supposed to be magnetars. But an 
alternative suggestion is that they are normal-field pulsar- 
like s t ars which are in an accretion propeller p hase iAlparl 
l200lt IChatteriee. Hernauist fc Naravanl l200of) . The very 
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difficulty in the later view point is to reproduce the 
irregular bursts, even with peak luminosity ~ 10 7 LEdd 
(SGR 0526-66; LEdd the Eddington luminosity). Both the 
possibilities of the bombardments of comet-like objects 
(e.g., strange planets, Xu 2005) to bare strange stars and of 
giant quakes in solid quark stars could remove the difficulty 
during super-bursts since the interaction in quark matter 
should be very strong. Based on the calculation in Fig|5| 
it is conjectured that SGRs/AXPs could be quark stars 
with masses in order of Mq, since more energy would be 
released during their quakes as bursts. Other pulsar-like 
stars (compact center objects and dim thermal "neutron 
stars"), which have small bolometric radii and low surface 
temperatures, could be quark stars with lower masses. 
Actually, there could be observations as well as theoretical 
arguments which do not favor the magnetar idea. (i). The 
superstrong field of magnetars are supposed to be created 
by MHD-dynamo action of rapid rotating protoneutron 
stars with spin period < 3 ms. The Poynting flux and 
the relativistic particle ejection of such a star should 
power effectively the supernova remnants. Such energetic 
remnants are e xpected in magnetar models, but had not 
been detected JVink fc Kuiperll2006l) . (ii). Dust emission 
around pulsar- like stars (e.g., AXPs) was proposed to test 
observationally the p ropeller sc e narios of quark stars with 
Spitzer or SCUBA iKul l2005t 120061) . Actually a recent 
discovery of mid-infrared emission from a cool disk around 
an isolated young X-ray pulsar, 4U 0142+61, is reported 
iWang. Chakrabartv fc Kaplan! [200611 although it is still a 
matter of debate whether significant propeller torque of 
fallback matter acts on the center star, (iii) The pressure 
should be very anisotropic in a relativistic degenerate neu- 
tron gas in equilibrium with a background of electrons and 
protons when the magnetic field is stronger than the critical 
field. The vanishing of the equatorial pressure of the gas 
would result in a transver se collapse, and a stable magneta r 
could then be unlikely iMartnez. Roias fc Cuestal [2003). 
(iv). It is still a matter of debate whether the absorption 
features in SGR 1806-20 can b e interpreted by proton o r 
electron cyclotron resonance lIXu. Wang fc Qiaol 12003^ . 
The field is only ~ 5 x 10 11 G in the c o ntext of electron 
cyclotron lines, (v). iMalov fc MachabeiH (120041) suggested 
that the persistent X-ray emission of AXPs/SGRs could 
be originated from the cyclotron mechanism acting near 
the surface of a star with normal field ~ 10 12 G, while 
short-time-scale cataclysmic events on the neutron star 
could lead to the bursts. Considering these criticisms, we 
think that alternative ideas for understanding AXPs/SGRs 
phenomena are very necessary. 

How to test this starquake model and the magnetar 
model for AXPs/SGRs? Because of the change of mass 
momentum, gravitational wave radiation may accompany 
a giant starquake of a solid quark star, but such ra- 
diation may not be significant in the magnetar model, 
at least the timescale and other feature of the waves 
should be very different. The observed QPOs could also 
hint the nature of superflares of SGRs. A detail inves- 
tigation on stellar torsional vibration and a comparison 
between the QPO-characters in starquake and magnetar 
models (e.g.. [Gl arnjjejiak^ I _Sam uelsson fc Anderssonll2006t 
iBeloborodov fc Thompsonl l2006) are very necessary. In our 
model for SGRs, glitches may occur during gravity-induced 



superflares. This could be applied to distinguish the models. 
Interestingly, an glitch of Av/v = 4.2 x 10 -6 was discovered 
in IE 2259+586, wh i ch pr eceded the the 2002 out-burst ac- 
tivity llWoods et al.ll2004l) . Actually, occasional 4 glitches 
had been detected in 3 AXP s; sometimes they were asso- 
ciated with radiative events (iKaspil I2006T) . However, only 
gravitational energy change is calculated in this paper. It 
is worth noting that stress energy developed in solid quark 
stars should be another kind of free energy. No significant 
frequency glitch occurs if a superflare originates from the re- 
lease of stress energy. How much is the stress energy during 
the evolution of a solid quark star? Is it comparable to the 
gravitational energy? These are questions to be answered in 
the future. 

A pulsar-like star could be monopole-charged electri- 
cally, due to either the global structure of cur rent flows 
in pulsar mag netospheres jXu^C ui fc Qiadl2006l) or maybe 
other reasons. IWeber et all J2OO6I) showed that, depending 
of the amount of electric charge, the star's structure and 
specifically the mass-radius relationship might be drastically 
modified. A loss of the electricity could also increase the pa- 
rameter e (i.e., the tangential pressure becomes higher and 
higher than the radial one as the electricity decreases). A 
starquake would occur as the star discharges, which we have 
not included in calculations. 

Astrophysical links between AXPs and SGRs? The per- 
sistent X-ray emission from AXPs/SGRs are alternatively 
suggested to be accretion-p owered (i.e., accretor with con- 
ventional magnetic fields, lAlparlEoOlh . But how to re- 
produce naturally SGR-like bursts in this scenario? Quark 
stars with low masses (<C Mq) are self-confined by the 
strong interactions between quarks, whereas gravitation- 
binding can not be negligible for quark stars with much 
higher masses (~ Mq). This results in an approximate re- 
lation of M oc R 3 for low masses but violation for higher 
masses in the mass(Af)-radius(_R) diagram. It is well known 
that the stellar radius decreases as the mass increases for 
pure gravitation-confined Fermion stars (e.g., white dwarfs 
with state of perfect electron gases). An accreting solid 
quark star may undergo blazing quakes if it has a mass of 
> Mq when gravitation-binding dominates, since, in this 
case, the gravitation-induced shear force becomes stronger 
and stron ger as the star accretes. Therefore, it is proposed 
JXul2006ft that SGRs/AXPs might be solid quark stars with 
~ (1 — 2)Mq, while other pulsar- like stars could be of low- 
mass. In order to make sense about the values of mass and 
radius of gravitation-dominated quark stars, we calculate 
these two parameters for fluid quark stars, using simply 
an equation of state of Eq.@. It is shown in Fig. 01 that 
the difference between the mass of quark stars with maxi- 
mum mass (AM/AR = 0) and that with maximum radius 
(AR/AM = 0) could be AM ~ 0.1M Q for possible bag con- 
stant (60 — 110) MeV-fm 3 . The corresponding radius differ- 
ence is about AR ~ 0.4 km. The maximum gravitational- 
energy release during successive quakes of an accreting solid 
quark star from the point of " AR/AM — 0" to the point of 
"AM/AR = 0" could typically be B quake ~ G{ M( ^i + / f^ I) - 

~ 10 52 erg. But the total energy release for the persis- 
tent X-ray emission is E pcl - sistcnt ~ GMAM/R ~ 10 erg. 
For AXPs/SGRs with persistent X-ray luminosity L x ~ 10 35 
erg/s, the lifetime of such sources could be E pels i EtC nt / L x ~ 
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emission to be fallback accretion-powered, for AXPs/SGRs. 
Those problems in the magnetar models could be circum- 
vented if AXPs/SGRs are accreting solid quark stars with 
masses ~ (1 — 2)Mq. The rigidity of solid quark stars should 
be high enough for strong quakes to occur, and the energy 
budget would not be a problem. 

Acknowledgments: RXX thanks Dr. Bing Zhang for a 
discussion of possible links between AXPs and SGRs. This 
work is supported by NSFC (10573002) and the Key Grant 
Project of Chinese Ministry of Education (305001). 
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Figure 4. An illustration of stellar mass (solid lines, in Mq) and 
radius (dashed lines, in km) of gravitation-dominated quark stars, 
as a function of the bag constant. The masses as well as the radii 
of quark stars with maximum radius (di?/dM = 0) or maximum 
mass (dM/dR = 0) are shown. It is evident that the mass and 
radius differences could be about O.IMq and 0.4 km, respectively. 



10 yrs. This means that it would cost nearly the Hubble 
time to increase the stellar mass to the maximum value (so 
that the star may collapse to a black hole) if the accretion 
rate is a constant. 

Radiative mechanisms of star-quake-induced flares. Be- 
cause of the starqua kes of neutron stars, self-induction elec- 
tric field is created llTho mpson. Lvutikov fc Kulkarnill2002t 
iBeloborodov fc Thompson! I2006T) . The strong electric field 
could initiate avalanches of pair creation in the magneto- 
sphere and certainly accelerate particles, resulting in high 
energy bursts observed. Similar mechanism would also work 
when starquakes occur in solid quark stars. The only differ- 
ence should be that no ions can be supplied from the sur- 
face of bare quark stars (i.e., lepton-dominated plasma forms 
above quark surfaces). The energy of oscillations excited by 
starquakes could be transported to the plasma corona by, 
e.g., Alfv en waves, being similar to the case of heating the so- 
lar wind IITull988l) . It is interesting to know if such Alfvenic 
fluctuations, originated from stellar torsional oscillations, 
could result in the observed QPOs. In the conventional mag- 
netar model, the field should be at least B ~ 8 x 10 1B G in 
order to reproduce the superflare of SRG 1806-20, if the ef- 
ficiency of transforming magnetic energy to the burst emis- 
sion during the flare is ~ 10~ 2 . This low limit means that 
magnetic field dominates (i.e., to control the motion by mag- 
netic field rather by crust) in the part of crust with density 
p < Pmax ~ 2.7 x 10 11 g/cm 3 . Note that p max is in a same 
order of the density of neutron drip pdrip = 4.3 x 10 11 g/cm 3 . 
It is possible (or not) that starquake could occur in the out- 
crust (with density p < pdrip) of normal neut ron stars, bu t 
the disadvantage is that the shear modulus l|Fuchsl [l92fi'). 
p. oc Z 2 , in the deep crust might not be high enough to crack, 
since the atomic charge Z becomes smaller and smaller in 
deeper crust due to neutronization. However, for solid quark 
stars, starquakes could certainly exist because of hight shear 
modulus as well as high density, even if the stars have field 
as strong as 10 16 G. 

In conclusion, starquakes of solid quark stars are pro- 
posed for the soft 7-ray superflares, with the persistent X-ray 
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